Emergence of Quantum Chaos in Quantum Computer Core and How to Manage It 
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We study the standard generic quantum computer model, which describes a realistic isolated quan- 
tum computer with fluctuations in individual qubit energies and residual short-range inter-qubit 
couplings. It is shown that in the limit where the fluctuations and couplings are small compared 
to one-qubit energy spacing the spectrum has a band structure and a renormalized Hamiltonian is 
obtained which describes the eigenstate properties inside one band. The studies are concentrated on 
the central band of the computer ( "core" ) with the highest density of states. We show that above a 
critical inter-qubit coupling strength, quantum chaos sets in, leading to quantum ergodicity of the 
computer eigenstates. In this regime the ideal qubit structure disappears, the eigenstates become 
complex and the operability of the computer is quickly destroyed. We confirm that the quantum 
chaos border decreases only linearly with the number of qubits n, although the spacing between 
multi-qubit states drops exponentially with n. The investigation of time-evolution in the quantum 
computer shows that in the quantum chaos regime, an ideal (noninteracting) state quickly disap- 
pears and exponentially many states become mixed after a short chaotic time scale for which the 
dependence on system parameters is determined. Below the quantum chaos border an ideal state 
can survive for long times and be used for computation. The results show that a broad parameter 
region does exist where the efficient operation of a quantum computer is possible. 

PACS numbers; 03.67.Lx, 05.45.Mt, 24.10.Cn 



I. INTRODUCTION 

During the last decade, a remarkable progress has been 
achieved in the fundamental understanding of the main 
elements necessary for the creation of a quantum com- 
puter. Indeed, as stressed by Feynman classical 
computers have tremendous problems to simulate very 
common quantum systems, since the computation time 
grows exponentially with the number of quantum parti- 
cles. Therefore for such problems it is natural to envision 
a computer composed from quantum elements (qubits) 
which operate according to the laws of quantum me- 
chanics. In any case, such devices will be in a sense un- 
avoidable since the technological progress leads to chips 
of smaller and smaller size which will eventually reach 
the quantum scale. At present a quantum computer is 
viewed as a system of n qubits (two-level quantum sys- 
tems) , with the possibility of switching on and off a cou- 
pling between them (see the detailed reviews in [||-§|). 
The operation of such computers is based on reversible 
unitary transformations in the Hilbert space whose di- 
mension Nh = 2" is exponentially large in n. It has been 
shown that all unitary operations can be realized with 
two-qubit transformations This makes necessary 

the existence of a coupling between qubits. Any quan- 
tum algorithm will be a sequence of such fundamental 
transformations, which form the basis of a new quantum 
logic. 

An important next step was the discovery of quantum 
algorithms which can make certain computations much 
faster than on a classical computer. The most impres- 
sive is the problem of factorization of large numbers in 



prime factors, for which Shor constructed a quantum 
algorithm which is exponentially faster than the classical 
ones. It was also shown by Grover ||^ that the searching 
of an item in a long list is parametrically much faster on 
a quantum computer. The recent development of error- 
correcting codes |^,|lO| showed that a certain amount of 
noise due to external coupling could be tolerable in the 
operation of a quantum computer. 

All these exciting developments motivated a great 
body of experimental proposals to effectively realize such 
a quantum computer. They include ion traps nu- 
clear magnetic resonance systems , nuclear spins with 
interaction controlled electronically or by laser 

pulses quantum dots jl^. Cooper pair boxes 
optical lattices [|8| and electrons floating on liquid he- 
lium ||l9|. As a result, a two-qubit gate has been ex- 
perimentally realized with cold ions |20| , and the Grover 
algorithm has been performed for three qubits made from 
nuclear spins in a molecule |2l|| . However, to have a quan- 
tum computer competitive with a classical one will re- 
quire a much larger number of qubits. For example, the 
minimal number of qubits for which Shor's algorithm will 
become useful is of the order of n = 1000 As a result, 
a great experimental effort is still needed on the way to 
quantum computer realization. 

A serious obstacle to the physical realization of such 
computers is the quantum decoherence due to the cou- 
plings with the external world which gives a finite life- 
time to the excited state of a given qubit. This ques- 
tion has been discussed by several groups for different 
experimental qubit realizations P,|6 p^ , |2^ . The effects 
of decoherence and laser pulse shape broadening were 
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numerically simulated in the context of Shor's algorithm 
p^ , p5t , and shown to be quite important for the operabil- 
ity of the computer. However, in a number of physical 
proposals, for example nuclear spins in two-dimensional 
semiconductor structures, the relaxation time due to this 
decoherence process can be many orders of magnitude 
larger than the time required for the gates operation 
P, p"3p4|j2^ , so that there are hopes to manage this ob- 
stacle. 

Here we will focus on a different obstacle to the 
physical realization of quantum computers that was 
not stressed up to now. This problem arises even if 
the decoherence time is infinite and the system is iso- 
lated/decoupled from the external world. Indeed, even in 
the absence of decoherence there are always imperfections 
in physical systems. Due to that the spacing between the 
two states of each qubit will fluctuate in some finite de- 
tuning interval S. Also, some residual static interaction J 
between qubits will be unavoidably present (we remind 
that an inter-qubit coupling is required to operate the 
gates). Extensive studies of many-body interacting sys- 
tems such as nuclei, complex atoms, quantum dots and 
quantum spin glasses p6|-p5t have shown that generically 
in such systems the interaction leads to quantum chaos 
characterized by ergodicity of the eigenstates and level 
spacing statistics as in Random Matrix Theory (RMT) 
1^,^. In a sense the interaction leads to dynamical 
thcrmalization without coupling to an external thermal 
bath. If the quantum computer were in such a regime, its 
operability would be effectively destroyed since the non- 
interacting multi-qubit states representing the quantum 
register states will be eliminated by quantum ergodicity. 

In this respect, it is important to stress that unavoid- 
ably the residual interaction J will be much larger than 
the energy spacing A„ between adjacent eigenstates of 
the quantum computer. Indeed the residual interaction 
J is relatively small so that all Nh computer eigenener- 
gies are distributed in an energy band of size AE ~ uAq, 
where Aq is the average energy distance between the two 
levels of one qubit and n is the total number of qubits 
in the computer. As a consequence, the spacing between 
multi-qubit states is A„ « AE/Nh ^ nAo2-" < Aq. 
Let us consider a realistic estimate for A„ and J for 
the case with n = 1000 as required for Shor's algo- 
rithm to be useful. For Aq ~ 1 K, which corresponds 
to the typical one-qubit spacing in the experimental pro- 
posals |l3| , [l4| , the multi-qubit spacing becomes A„ ~ 
10^ X 2-10 Ao - 10-298 K. This value will definitely be 
much smaller than any physical residual interaction. In 
the case of the proposal ||l^, for example, with a dis- 
tance between donors of r = 200 A and an effective Bohr 
radius of as = 30 A ( Eq.(2) of Q), the couphng be- 
tween qubits (spin-spin interaction) is J ~ Aq ~ 1 K. 
By changing the electrostatic gate potential, the effective 
electron mass can be modified up to a factor of two. Since 



J oc {r/aB)^/"^ exp{—2r/aB)/a,B: and is inversely pro- 
portional to the effective mass, this gives a minimal resid- 
ual spin-spin interaction of J ^ 10-^ K ^ A„. In this 
situation, one would naturally/naively expect that level 
mixing, quantum ergodicity of eigenstates and chaos are 
unavoidable since the interaction is much bigger than the 
energy spacing between adjacent levels ( J ^ A„). 

In spite of this natural expectation, it was shown re- 
cently in that in the quantum computer the quantum 
chaos sets in only for couplings J exponentially stronger 
than A„. In fact, it was shown that the critical coupling 
Jc for the transition to quantum chaos decreases only 
linearly with the number of qubits n (for short-range 
inter-qubit coupling): Jc ~ Ag/n. This result opens 
a broad parameter region where a quantum computer 
can be operated below the quantum chaos border, when 
noninteracting multi-qubit states are very close to the 
exact quantum computer eigenstates. For example, at 
n = 1000 and Aq '--^ 1 K, the critical residual interaction 
is Jc '--^ 1 mK, compatible with the proposal discussed 
above II . 

In the present paper, we study in more details the tran- 
sition to chaos and how it affects the time evolution of the 
system. The effects of residual interaction in the pres- 
ence or absence of fine fluctuations of individual qubit 
energy spacing are analyzed in great detail. The paper 
is composed as follows. In the next section we describe 
the standard generic quantum computer (SGQC) model, 
introduced in |38|. In section III, we present the result 
of numerical and analytical studies of eigenenergies and 
eigenstate properties of this model. Section IV is devoted 
to the analysis of the time evolution of this system, and 
the typical time scales for the development of quantum 
chaos are presented as a function of the system param- 
eters. We end by some concluding remarks in the last 
section. 

II. STANDARD GENERIC QUANTUM 
COMPUTER MODEL 

In p8| the standard generic quantum computer 
(SGQC) model was introduced to describe a system of n 
qubits containing imperfections which generate a resid- 
ual inter-qubit coupling and fluctuations in the energy 
spacings between the two states of one qubit. The Hamil- 
tonian of this model reads: 

H = J2T,a^ + Y.J^,<7^<7^, (1) 

i i<j 

where the ai are the Pauli matrices for the qubit i and 
the second sum runs over nearest-neighbor qubit pairs on 
a two-dimensional lattice with periodic boundary condi- 
tions applied. The energy spacing between the two states 
of a qubit is represented by F.^ randomly and uniformly 
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distributed in the interval [Aq - 5/2, Aq + 5/2]. The de- 
tuning parameter 5 gives the width of the distribution 
near the average value Aq and may vary from to Aq. 
Fluctuations in the values of Ti appear generally as a re- 
sult of imperfections. For example, in the frame of the 
experimental proposals |13 1^, the detuning 5 will ap- 
pear for nuclear spin levels as a result of local magnetic 
fields and density fluctuations. For electrons floating on 
liquid helium it will appear due to fluctuations of 
the electric field near the surface. The couplings Jij 
represent the residual static interaction between qubits 
which is always present for reasons explained in the intro- 
duction. They can originate from spin-exciton exchange 
[^,n. Coulomb interaction dipole-dipole interac- 
tion [Q, etc... To catch the general features of the dif- 
ferent proposals, we chose Jij randomly and uniformly 
distributed in the interval [—J, J]. We note that a sim- 
ilar Hamiltonian, but without coupling/detuning fluctu- 
ations, was discussed for a quantum computer based on 
optical lattices This SGQC model describes the 

quantum computer hardware, while the gate operation 
in time should include additional time-dependent terms 
in the Hamiltonian (|l|) and will be studied separately. At 
J = the noninteracting eigenstates of the SGQC model 
can be presented as >= \ai, an > where = 0, 1 
marks the polarization of each individual qubit. These 
are the ideal eigenstates of a quantum computer, and 
we will call them quantum register states. For J ^ 0, 
these states are no longer eigenstates of the Hamiltonian, 
and the new eigenstates are now linear combinations of 
different quantum register states. We will use the term 
multi-qubit states to denote the eigenstates of the SGQC 
model with interaction but also for the case J = 0. 
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FIG. 1. Density of multi-qubit states of (|l|) as a function of 
total system energy E for J = 0. Here n — 16 and 5/Ao — 0.2. 
The two extreme bands at E/ Ao ~ ±16 contain only one state 
and are not seen at this scale. 

While in |38| the main studies were concentrated on 



the case where S is relatively large and comparable to 
Aq, here we will focus on the case 6 <C Aq, which cor- 
responds to the situation where fluctuations induced by 
imperfections are relatively weak. In this case, the unper- 
turbed energy spectrum of ([^) (corresponding to J = 0) 
is composed of n -I- 1 well separated bands, with inter- 
band spacing 2Ao. An example of the density of multi- 
qubit states pn = 1/A„ in this situation is presented in 
Fig.|l|. Since the Ti randomly fluctuate in an interval of 
size S, each band at J = 0, except the extreme ones, 
have a Gaussian shape with width y/nS. The number 
of states inside a band is approximately Nn/n, so that 
the energy spacing between adjacent multi-qubit states 
inside one band is exponentially small (5„ ~ n'^/^2~"5), 
in line with the general estimate in Section I. 

In the presence of a residual interaction J ~ 5, the 
spectrum will still have the above band structure with 
exponentially large density of states. For J ~ (5 <C Aq, 
the interband coupling is very weak and can be neglected. 
In this situation, the SGQC Hamiltonian (|^) is to a good 
approximation described by the renormalized Hamilto- 
nian Hp = Yi^t\PkH Pk where i\ is the projector on 
the fc*'' band, so that qubits are coupled only inside one 
band. We will thereafter concentrate our studies on the 
band nearest to i? = 0. For an even n this band is cen- 
tered exactly at i? = 0, while for odd n there are two 
bands centered at _E = ±Ao, and we will use the one 
at E = — Aq. Such a band corresponds to the highest 
density of states, and in a sense represents the quantum 
computer core. It is clear that quantum chaos and ergod- 
icity will first appear in this band, which will therefore 
set the limit for operability of the quantum computer. 
Inside this band, the system is described by a renormal- 
ized Hamiltonian Hp which depends only on the number 
of qubits n and the dimensionless coupling J /5. 



III. QUANTUM COMPUTER EIGENENERGIES 
AND EIGENSTATES 

The first investigations in showed that the quan- 
tum chaos border in the SGQC model (|^) corresponds to 
a critical interaction Jc given by: 



Jr. 



C6 

n 



(2) 



where C is a numerical constant. This border is expo- 
nentially larger than the energy spacing between adjacent 
multi-qubit states A„. The physical origin of this differ- 
ence is due to the fact that the interaction is of a two- 
body nature. As a result, one noninteracting multi-qubit 
state \ipi > has nonzero coupling matrix elements only 
with 2n other multi-qubit states. In the basis of quantum 
register states {ipi >^ the Hamiltonian is represented by a 
very sparse nondiagonal matrix with only 2n -I- 1 nonzero 
matrix elements by line of length Nh = 2". For 5 Aq 
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all these transitions take place in an energy interval B 
of width of order 6Ao. Therefore the energy spacing be- 
tween directly coupled states is Ac ~ B/2n w 3Ao/n. 
According to the studies of quantum chaos in many-body 
systems p7| , |30| - |35| , ^ , the transition to chaos takes place 
when the matrix elements become larger than the en- 
ergy spacing between directly coupled states. This gives 
J > Ac which leads to the relation (||). For the case 
6 <C Aq on which we focus here, still in the renormalized 
Hamiltonian Hp the number of nonzero matrix elements 
in one line is of the order of n, and B S, so that 
Ac^ 5/n, that leads to the result (0) Q. 




1 2 33"* 

FIG. 2. Transition from Poisson to WD statistics for the 
renormalized Hamiltonian of SGQC model in the central 
band. The statistics is obtained for the states in the mid- 
dle of the energy band (±6.25% around the center) for 
n=16 : J/S = 0.05,7) — 0.99 (dashed line histogram); 
J/S — 0.32, 77 = 0.047 (full line histogram). Full curves show 
the Poisson distribution Pp{s) and the Wigner-Dyson distri- 
bution Pw{s); Nd = 8, iVs > 1.2 x 10*. 

The transition to quantum chaos and ergodicity can be 
clearly seen in the change of the spectral statistics of the 
system. One of the most convenient is the level spacing 
statistics P{s), which gives the probability to find two ad- 
jacent levels whose spacing is in [s, s + ds]. Here s is the 
energy spacing measured in units of average level spac- 
ing. It is well known that while the average density of 
states is not sensitive to the presence or absence of chaos, 
the fluctuations of the energy spacings between adjacent 
levels around the mean value, determined by -P(s), are 
sensitive to it. In the presence of chaos, eigenstates are 
ergodic, overlap of wavefunctions gives a finite coupling 
matrix element between nearby states and the spectral 
statistics P{s) follows the Wigner-Dyson (WD) distri- 
bution Pw{s) — (7rs/2) exp(— 7rs^/4) typical for random 
matrices. This distribution Pwis) shows level repulsion 
at small s, due to the fact that overlap matrix elements 
between adjacent levels tend to move them away from 
each other. On the contrary, in the integrable case at 



J <C Jc, the overlap coupling matrix element between 
nonergodic states is very small. As a result, energy levels 
are uncorrelated and P(s) follows the Poisson distribu- 
tion Pp{s) = exp(— s) known to be valid for integrable 
one-particle systems 

In the SGQC model, we expect a transition from Pp{s) 
at small J to Pw{s) above the quantum chaos border 
(^. An example of such a transition is shown in Fig. 2. 
To decrease the statistical fluctuations we averaged over 
several independent realizations of the F,; and in (|^), 
which is the standard procedure used in Random Matrix 
Theory | p6| , p7t . We used up to No ~ 5x 10'' realizations 
so that the total statistics 1.5 x 10^ > Ns > 1.2 x 10''. 
It is interesting to note that in the limit J/S — > 00 ((5 <C 
J <C Ao) the system remains in the regime of quantum 
chaos with WD statistics Q], as is illustrated in Fig. 3. 
This means that in the absence of individual qubit energy 
fluctuations, the residual coupling alone leads to chaotic 
eigenstates. 
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FIG. 3. Level spacing statistics for the renormalized Hamil- 
tonian of SGQC model in the central band for S — 0. The 
statistics is obtained for the states in the middle of the en- 
ergy band (±6.25% around the center) for n=15 : 77 = 0.023 
(histogram). Full curves show Pp{s) and Pw{s); Nd = 20, 
Ns > 1.6 x lO''. 

To characterize the variation of P{s) from one limiting 
distribution to another it is convenient to use the parame- 
ter r, = J^" (P(s) - Pw {s))ds/ J^" {Pp (s) - Pw {s))ds ^ , 
where sq = 0.4729... is the intersection point of Pp{s) 
and Pwis). In this way Pp{s) corresponds to 77 = 1, and 
Pw{s) to 1] — 0. The studies of different systems has 
already shown that this parameter characterizes well the 
transition from one statistics to the other |31 3^,^,^. 
Indeed, according to the data of Fig. 4, 77 changes from 
1 at small J to 77 « at large J. To characterize this 
transition, we chose the critical value Jc by the condition 
ry(Jc) = 0.3. The dependence of 77 on the rescaled cou- 
pling strength J/ Jc shows that the transition becomes 
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sharper and sharper when n increases (Fig. 4). 




FIG. 4. Dependence of rj on the rescaled coupling strength 
J/Jc for the states in the middle of the energy band for 
n = 6(*), 9(0), 12 (triangles), 15(squares), 16 (diamonds). 

The dependence of the critical coupUng strength Jc on 
the number of qubits n is shown on Fig. 5. It clearly shows 
that this critical strength decreases linearly with n and 
follows the theoretical border (^ with C « 3. For com- 
parison on the same figure we also show the dependence 
of the multi-qubit spacing A„ (computed numerically) 
on n. It definitely demonstrates that Jc 3> A„. 
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FIG. 5. Dependence of \og{Jc/5) (diamonds) and 
log( Jcs/5) (triangles) versus log(n); the variation of the scaled 
multi-qubit spacing (log(An/(S)) with log(n) is shown for 
comparison (-I-). Dashed line gives the theoretical formula 
Jc ~ C5/n with C = 3.3; the solid line is Jcs — OAlS/n; the 
dotted curve is drawn to guide the eye for {+). 

The transition in the level spacing statistics reflects 
a qualitative change in the structure of the eigenstates. 
While for J <C Jc the eigenstates are expected to be very 
close to the quantum register states jV'i >, for J > Jc 



each eigenstate \(f>rn > becomes a superposition of an ex- 
ponential number of states \tjji >. It is convenient to 
characterize the complexity of an eigenstate \(j)m > by the 
quantum eigenstate entropy Sq = — log2 W^m, 

where Wim is the quantum probability to find the quan- 
tum register state {tpi > in the eigenstate {(pm > of the 
Hamiltonian {Wim — \ < ipi\4'm > P)- In this way Sq = 
if \(j>rn > is one quantum register state ( J = 0), iSg = 1 if 
> is equally composed of two >, and the maxi- 
mal value is Sq — n if all 2" states contribute equally to 
\4>m >■ We average Sq over the states in the center of 
the energy band and Njj realizations of F^ and Jij . 
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FIG. 6. Dependence of the quantum eigenstate entropy Sq 
on J/Jc for n — 6(*), 9 (o), 12 (triangles), 15 (squares), 16 
(diamonds); 1.5 x 10" > Ns > 1.2 x lO''. 
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FIG. 7. Same as Fig. 6 but on a larger scale. 



The variation of this average Sq as a function of J for 
different values of n is shown on Figs. 6, 7. It shows that 
indeed the entropy Sq grows with J until it saturates to 
a large value corresponding to an exponential number of 
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mixed states. These data show that the critical couphng 
Jcs at which Sq = 1 (two states mixed) is proportional to 
Jc- Indeed, Fig. 7 shows a small dispersion near 5"^ = 1 
when n changes from 6 to 16, while A„ varies by three 
orders of magnitude. This is confirmed by the data on 
Fig. 5, which give Jcs ~ O.lSJc ~ QA5/n. This result 
is in agreement with the results |38[ obtained by direct 
diagonalization of the SGQC model at 5 ^ Ao (lower 
insert in Fig. 2 of ||||). 

The quantum eigenstate entropy Sq characterizes the 
global properties of the eigenstates, while a more detailed 
information about them can be obtained from the local 
density of states pw introduced by Wigner fi^ : 



Pw {E — Ei 



W,r,^5{E - E„ 



(3) 



iog(r/8) 




-0.5 0.0 

log(J/ 8) 

FIG. 8. Dependence of the Breit-Wigner width F on the 
coupling strength J for n — 15 for the states in the middle 
of the energy band. The straight lines show the theoretical 
dependence 1^ with F = 1.3J^n/S and the strong coupling 
regime with F ~ J; Nn — 20. Lower insert: example of 
the local density of states pw (^ for J/5 — 0.08; the full line 
shows the best fit of the Breit-Wigner form (^ with F = 0.105. 
Upper insert: example of the local density of states pw ^ for 
J/S = 0.4; the full line shows the best Gaussian fit of width 
F = 0.645. 

The function pw characterizes the average probability 
distribution of Wim (see a numerical example in Fig. 3 
of Q). For moderate couphng strength, pw is well 
described by the well-known Breit-Wigner distribution 
Pw = Pbw- 



Pbw{E — Ei) = 



(4) 



where F is the width of the distribution. This expression 
is valid when F is smaller than the bandwidth (F < y/nS) 
and many levels are contained inside this width. In this 
regime, the Breit-Wigner width F is given by the Fermi 



golden rule: F = 2TrU^pc, where Us is the root mean 
square of the transition matrix element and pc is the 
density of directly coupled states. The validity of this 
formula was well checked in many-body systems with 
quantum chaos | p8| , ^ ,p4 37 1. In our case Us ^ J and 
Pc ~ n/6, so that: 



F - 



(5) 



This dependence is confirmed by the data on Fig. 8. 
However, for large J, when F > ^/n5^ the shape of pvi' 
becomes non-Lorentzian and is well fitted by a Gaus- 
sian distribution. The width of this modified distribu- 
tion grows like F ~ J. This scaling naturally appears in 
the limit (5 = 0, J <C Aq, since the noninteracting part 
of the Hamiltonian is simply a constant commuting with 
the perturbation. The change from one dependence to 
the other takes place for J > S/v}^^. Above this limit F 
is still weakly dependent on the number of qubits n. We 
expect that for J ^ S the energy width of one band is 
F ~ J\/n (effective frequency of n Rabi frequencies with 
random signs), and have checked numerically this law for 
(5 = (data not shown). 




FIG. 9. Melting of the quantum computer core generated 
by the inter-qubit coupling. Color represents the level of 
quantum eigenstate entropy Sq, from bright red {Sq ~ 12) 
to blue {Sq — 0). Horizontal axis is the scaled energy E/S 
of the computer eigenstates in the central band counted from 
the band bottom to the top {E/S ~ it^'n). Vertical axis is the 
value of J/5, varying from to 0.5. Here n = 16, Jc/S = 0.22, 
and one random realization is chosen. 

According to the results obtained from many-body sys- 
tems |]33| , the number of quantum register states mixed 
inside the width F is of the order of Fp„, and is ex- 
ponentially large. This however assumes that J > Jc 
and the system is already in the quantum chaos regime. 
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In this case the quantum eigenstate entropy Sq is large 
{Sq « log2(rp„) ~ n) and the operabihty of the com- 
puter is quickly destroyed, since many quantum register 
states become mixed. The pictorial view of the quantum 
computer melting is shown on Fig. 9. This image is quali- 
tatively similar to the one in (Fig. 5 there), which was 
obtained for the SGQC model at ^ = Aq. In Fig. 9 the 
melting goes in a smoother way since all the states belong 
to the same central band (quantum computer core). 




FIG. 10. Quantum chaos in the quantum register: color 
represents the value of the projection probability Wim of 
the quantum register states on the eigenstates of the Hamil- 
tonian, from bright red (maximal vaIue)to blue (minimal 
value). Horizontal axis corresponds to 150 quantum register 
states, the vertical axis represents the 150 computer eigen- 
states (both ordered in energy). Here n — 16, J/5 = 0.4 
{J/5 > Jc/5 — 0.22), and one random realization is chosen. 




FIG. 11. Same as Fig. 10 below the quantum chaos border, 
J/5 = 0.001 (J/5 < Jcs/5 = 0.026). 




FIG. 12. Same as Fig.lO for J/5 = 0.01 
{J/5 ~ Jcs/5 = 0.026) 

The effect of quantum chaos melting in the quantum 
register representation is shown on Fig.lO for J > Jc- 
The ideal register structure is manifestly washed out. 
On the contrary, below the chaos border (J < Jc), only 
few quantum register states are mixed. For comparison, 
Fig. 11 shows the same part of the register in the regime 
J ^ Jcs (no mixing of states) and Fig. 12 in the regime 
J ^ Jcs (few states are mixed) . 

IV. TIME EVOLUTION IN THE SGQC MODEL 

In the previous section we determined the properties 
of eigenstates of the quantum computer in the presence 
of residual inter-qubit coupling. In the presence of this 
coupling the quantum register states > are not any 
more stationary states, and therefore it is natural to an- 
alyze how they evolve in time. Indeed, if at time < = 
an initial state is \x{t = 0) >= |i/'io > corresponding to 
the quantum register state io, then with time the prob- 
ability will spread over the register and at a time t the 
projection probability on the register state \4>i > will be: 

Fuoit) ^ \ < ^Mt) > \' If.. 

= E™,™' A^n.AlmA*„^'A^,rn' exp(i(i?„, - E^)t), 

where Aim =< '4^i\4>m > and Em is the energy of the 
stationary state \(j>m > and we chose h = 1. For J <^ 
the probability Eigig{t) is very close to one for all times 
since the states are not mixed by the interaction. This 
means that all quantum register states {ipi > remain well 
defined, and the computer can operate properly. For J ^ 
Jcs, only few states > are mixed by the interaction, 
and Eigig {t) oscillates in time regularly around an average 
value of order 1/2. These oscillations are similar to the 
Rabi oscillations between two levels with frequency O ~ 
J. An example is presented in Fig. 13. In this regime, we 
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expect that error-correcting codes may efficiently 

correct tlie spreading over few quantum register states. 



0.6 - 




mined in the previous section. This estimate is very nat- 
ural in the Fermi golden rule regime, with Breit-Wigner 
local density of states (^), since Figig{t) is essentially the 
Fourier transform of the local density of states pw, and 
therefore decreases as exp{—Tt). We note that the decay 
in this regime was recently discussed in Q. According 
to our data, when F becomes comparable to the energy 
bandwidth y/nS, pw is close to a Gaussian distribution 
of width F, and its Fourier transform Figig{t) is also a 
Gaussian of width 1/F. Therefore in both regimes we 
expect the time scale t-^ for the decay of Fi^i^{t) to be 
^ 1/F. The data shown on Fig. 14 correspond to the 
saturation regime for large values of n, and the insert 
shows that ~ 1/F is still valid. In fact the curve for 
n = 16 in Fig. 14 is already close to the limiting decay 
curve &i d — Q (data not shown). 




FIG. 13. Time-dependence of the probability to re- 
main in the same quantum register state for n — 16, 
J = 0.01 ~ Jcs = 0.026 ( Jc/5 = 0.22) and one random 
realization is chosen. 





FIG. 15. Time-dependence of the quantum entropy S{t) for 
J/5 — 0.4 2> Jc/5; symbols are as in Fig. 14. Average is made 
over 200 initial states randomly chosen in the central band. 
Insert shows the same curves normalized to their maximal 
value. 



0.0 ' ■ ' ■ — ^y=#=^g-^gg~ijrtrT^6 

1 2 ^ t5 

FIG. 14. Time-dependence of the probability to remain in 
the same quantum register state for J/5 — 0.4 3> Jc/5. Data 
are shown for n = 16 (diamonds, Jc/5 = 0.22); n — 15 
(squares, Jc/5 = 0.24); n = 12 (triangles, Jc/5 = 0.28); 
n = 9 (circles, Jc/5 — 0.35); n — 6 (stars, Jc/5 — 0.59). 
Average is made over 200 states randomly chosen in the cen- 
tral band. Insert shows the chaotic time scale (defined 
by Figif^{T^) — 1/2) as a function of l/F; the straight line is 

= 1.27/F. 

For J > Jcs, quantum chaos sets in, and with time 
the probability spreads over more and more quantum 
register states until a quasi-stationary regime is reached 
where an exponentially large number of states are mixed. 
The probability Figig{t) drops approximately to zero, as 
shown on Fig. 14. The chaotic time scale for this decay 
can be estimated as ~ 1/F where F is the width deter- 



At the same time scale the quantum entropy 
S{t) is large but is still growing. It reaches its max- 
imal value on a larger time scale which seems inde- 
pendent of n. At this stage, an initial quantum regis- 
ter state is now spread over most of the register (Here 
S{t) = -I]i-fiio(^)log2^"oW)- This process is shown 
on Fig. 15. This maximal value of S{t) is approximately 
given by Sq (see Fig. 6) and accordingly decreases with 
decreasing J as is illustrated in Fig. 16. 

Fig. 17 illustrates this mixing process in the quantum 
register representation, evolving in time. The quantum 
computer hardware becomes quickly destroyed due to the 
inter-qubit coupling. It is necessary to decrease the cou- 
pling strength below the quantum chaos border to get 
well-defined quantum register states for t > 0, as is illus- 
trated in Fig. 18. 
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FIG. 16. Time-dependence of the quantum entropy S{t) 
for different values of J, n = 16 Jc/5 = 0.22, = 0.026, 

and one random realization is chosen: J/S = 0.001 <^ Jcs/S 
(disks); J/S = 0.01 < Jcs/S (crosses); J/S = 0.03 ^ Jcs/S 
(squares); J/S = 0.2 ^ Jc/S (diamonds); J/S = 0.4 > Jc/S 
(triangles) . Insert gives the probability to remain in the same 
quantum register state for the same values of J/S. Averages 
are made over 200 states randomly chosen in the central band. 




FIG. 17. Time explosion of quantum chaos in the quantum 
register: color represents the value of the projection prob- 
ability I < ipilxii) > P of an initial state on the quantum 
register states ordered in energy, from bright red (maximal 
value)to blue (minimal value). Horizontal axis corresponds 
to 150 states, the vertical axis to 150 time steps, from tS = 
to tS — 2. At tS — 0, the chosen initial state is the superposi- 
tion of two quantum register states. Here n — 16, J/S = 0.4 
{J/S > Jc/5 — 0.22), and one random realization is chosen. 




FIG. 18. Same as Fig. 17 below the quantum chaos border, 
J/S = 0.001 {J/S < Jcs/S = 0.026). 

The obtained data clearly show that exponentially 
many quantum register states become mixed after the 
finite chaotic time scale « l/F. 

V. CONCLUSIONS 

The results presented in this paper show that residual 
inter-qubit coupling can lead to quantum chaos and er- 
godic very complicated eigenstates of the quantum com- 
puter. We have shown that in this regime the quantum 
register states disintegrate quickly in time over an expo- 
nentially large number of states and the computer op- 
erability is destroyed. We determined the dependence 
of the chaotic time scale of this process on coupling 
strength J, detuning fluctuations S of one-qubit energy 
spacing, and number of qubits n. After this time 
the quantum computer hardware is melt. To prevent 
this melting one needs to introduce an efficient error- 
correcting code which operates on a time scale much 
shorter than and suppresses the development of quan- 
tum chaos. To avoid the quantum chaos regime dan- 
gerous for quantum computing, one should engineer the 
quantum computer in the integrable regime below the 
quantum chaos border Jc ~ 3d/n. It is important to note 
that this border decreases with the detuning 6, show- 
ing that imperfections do not all conspire against the 
operability of the computer. We stress again that the 
transition to quantum chaos is an internal process which 
happens in a perfectly isolated system with no coupling 
to external world. Nevertheless, since a decoherence can 
be viewed as a result of internal interactions in a larger 
system, the results presented here may also apply to this 
problem. 
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Our main conclusion is that although in the quantum 
chaos regime a quantum computer cannot operate for 
long, fortunately the border for this process happens to 
be exponentially larger than the spacing between adja- 
cent computer eigenstates, and therefore a broad param- 
eter region remains available for realization of a quantum 
computer. Another possibility is to operate the quantum 
computer in the regime of quantum chaos. However, here 
one should keep in mind that after the chaotic time scale 
Ty. the computer hardware will melt due to inter-qubit 
coupling and quantum chaos. Therefor, the computer 
operability in this regime is possible only if many gate 
operations can be realized during the finite time (in 
a sense it becomes similar to the decoherence time). It 
is clear that the most preferable regime corresponds to 
quantum computer operation below the quantum chaos 
border. 
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